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1. Introduction 

Recent mathematical worlc ([|l],0; an introduction for physicists is in 0) has produced 
constructions of the 'master field' for general large A^ gauge and matrix field theories. The 
essential difficulty of the large number of degrees of freedom in higher dimensional large A^ 
theories is dealt with by finding master fields which live in 'large' operator algebras such 
as the type IIi factor associated with a free group. Singer has advocated these ideas as 
the appropriate framework for treating the large A^ limit of general field theories such as 
four-dimensional QCD. [|] 

So far the explicit constructions have been for very simple solvable theories such as 
decoupled multi-matrix integrals or two-dimensional Yang-Mills theory. The constructions 
required as input correlation functions computed using other techniques (such as summing 
diagrams or the existing saddle-point methods) and thus are not explicit for physically 
more interesting theories. One way to treat such theories would be to find equations 
satisfied by master fields and solve them. This approach was worked on extensively in the 
early 80's (some references are @-^) and several such equations were proposed, identical 
to or derived directly from the original classical equations of motion. However, the lack of 
analytic techniques for working with master fields stopped progress in this direction. 

In this paper we will apply the new ideas in the large A^ stochastic approach of 
Greensite and Halpern. 0] They proposed interpreting the Langevin equation of stochastic 
quantization, 

^</>(r,x) = -^^^;^^[./>]+r7(r,x) (1.1) 

as an equation for a master field (p. The main advantage of large A^ in this approach is 
that one eliminates the second step of stochastic quantization, averaging over the noise ry, 
by using a 'master noise source,' a master field r] reproducing Gaussian correlations. The 
claim is then that the proper time evolution has a large r fixed point which is the master 
field, and Greensite and Halpern showed this in perturbation theory. 

Our main result is to show how this can be made precise using a precise definition 
of master field, and to show its equivalence to the standard factorized Schwinger-Dyson 
equations. A complete treatment is given only for a solvable case (the one-matrix model) 
but the ideas are general. Techniques for solving the equations we discuss in multi-matrix 
or higher dimensional models do not exist, but at least one can describe what one is looking 
for. 

Let us make two remarks on similar results which should also exist. First, the canonical 
formalism was also considered P,|^,|[ and it appears to us that an analogous but even 
simpler statement might be made - that the field (p{x) and its canonical conjugate H(x) 
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can be reproduced by master fields which (in some sense) satisfy the classical equation 
of motion (p = {H, (p} with the original Poisson bracket. The simplicity of this is very 
attractive but a major unsolved problem in this approach is to characterize which points 
in phase space (</>, H) are possible master fields (for example, ^ = 11 = clearly is not). 
Second, many workers and in particular Greensite and Halpern found it useful to eliminate 
the space-time dependence of the fields in favor of structure in the internal variables by 
the procedures of 'reduction' and 'quenching.' We have little to say about this; it is not 
even clear to us whether it is a step forward or backward. 

In section 2 we discuss master fields, in section 3 we review stochastic quantization, 
and in section 4 give our main results. In this version, we have commented on the relation 



to earlier work of Haan |]T^ and of Cvitanovic, Lauwers and Scharbach |jT^ in section 5 



and in section 6 state conclusions. 



2. Master fields 



We will consider a matrix model with fields Ma with a G 5" a set, and action S'[M]. 
Most of what we say applies for a field theory if we take S to include space-time, 6^,^ to 
include a Dirac delta function, and so forth. 

It is worth making a few general comments about master fields. We can regard the 
set of all vacuum expectation values 

(Tr [Ma.Ma, ...))= lim li [ DM e"^^!^^] Tr M^.M^, . . . (2.1) 

as a linear functional W{a) on the algebra of words generated by formal variables TOq, 
with a & S. The functional depends on S and if we consider more than one action, we 
would indicate this dependence as Ws{a). In the large A^ limit, no kinematic (independent 
of S) equality relations between expectations of different words are known (the finite A^ 
Mandelstam relations certainly have no obvious limit) and we claim that none exist. (See 
1^ for arguments in this direction, which we will pursue elsewhere.) Thus to represent the 
general functional W, we need to consider variables ttIq, satisfying no relations: the algebra 
is the free algebra on n = \S\ (cardinality of S) generators; call it JF„. 

The problem of finding a master field is the problem of finding operators M^ and 
trace Tr such that 

!> Ma.Ma^ ■■■Ma,= W{ma,ma, . . . m^J. (2.2) 
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It turns out that the problem cannot be solved using the usual Tr M = ^^ Ma and one 
must instead consider more general traces satisfying linearity and Tr AB = Tr BA (which 
implies Tr U~^AU = Tr A for a unitary U). 

The solution will only be unique up to unitary equivalence. For a one matrix integral, 
the following solution is a perfectly good one: M = X and Tr O = f dXp{X)0. We might 
consider others in which M was not diagonal. For example, the method of orthogonal 
polynomials provides a representation M = ^ Rn{en,n+i + en+i,n) and Tr O = J2i=i ^a- 
Since the inner product does not enter ( |2.2| ), we might even consider general similarity 
transformations SMiS~^ with S any invertible operator. 

The first criterion we must satisfy is that the Mq. form a representation of the algebra 
JF„. At the formal level we work at in this paper, J^^ has no structure, and this gives no 
constraint on the Ma- As we will see shortly however the Mq, are operators on a very large 
Hilbert space and it is easy to write formal expressions which have convergence problems 
and probably do not make sense. To write the equations satisfied by physical master fields, 
we will want to use concepts like M^ or dM/dx and it would be very useful to know when 
these do make sense. For typical (regulated) actions the operators Mq, are bounded and 
thus the framework in which to make precise definitions is the theory of operator algebras. 
The algebra Tn is a subalgebra of a C* algebra, the group algebra of the free group on n 
generators, F„. 

This theory provides a construction of the master field given the loop functional Ws 
for any field theory: the GNS representation. This construction puts all the information 
from W in the inner product and it is not so explicit, but it does demonstrate the existence 
of the master field as well as yielding the important result that not every loop functional 
W corresponds to a master field: it must satisfy positivity constraints. 

Another representation has been provided by the work of Voiculescu. [Q The simplest 
case and a very relevant one for stochastic quantization is the Gaussian master field. Let 
us take (MqM^) = Sap- Since we can make linear transformations M^ = R"Ma, this 
is general. The variables Mq are free random variables in the sense of Voiculescu: the 
expectation of any word is determined by the expectations of the individual matrices in 
a way analogous to the construction of a joint distribution of independent commuting 
random variables. 

These ideas lead to an explicit operator construction of the master field: 

Mq = ttQ + a; (2.3) 

where Cq and a* are operators defined by the following matrix elements. We define the 
'free Fock space on 5" to be a Hilbert space with orthonormal basis vectors labelled by an 



ordered list of zero or more elements aj G S. 

Q'Q.|Ctn5 Oin—li • • • 1 Oil) |Q^5 Q^nj Q^n— 15 • • • 1 Q^l/ 

aa|Q;n,an-i,- ■■,ai) = 5a,a„|an-l,- • •,«!) (2.4) 

ftal) = 0. 

They are similar to bosonic creation and annihilation operators but with two differences. 
First, they are free - the product of two operators associated with a and (3 ^ a satisfies 
no relations. Second, the usual symmetry factor for bosonic harmonic oscillators is absent. 
Thus, we have not [a, a*] = 1 but instead 

5^K,a*] = |^|-l + |)(|. (2.5) 

a 

The trace is not the standard one (which does not make sense here) but is defined by 

ir A={\A\). (2.6) 

It is a trace in that it satisfies the axioms for a trace, for example tr [A, B] = 0, as we will 
see. 



The simplest proof of ( p.2|) is diagrammatic. If we expand the product of the 
Mon = Oq, . + a* in 2'^ terms, we can associate a planar diagram with each non-zero term. 
Consider the example in the figure of (tr MbM^M^M^). We start with Mctf,=a\) which 
makes the state \a), and we interpret this as adding a line with the label a to the diagram. 
With Mq,j._^, there are in general two possibilities: we can always act with a* _ to create 
a new line labeled ctfc-i, and if a/c-i = ctfc we can annihilate the existing line with aQ,^_^ as 
well. The 'free' nature of the Fock space precisely reproduces the planarity constraint on 
the diagrams. After applying Mct-^^, we must be left with no lines to have a matrix element 
with (|. This proves ( |2.2| ), and since Tr is a trace, so is tr . 

This is a special case of the 'free product of distributions,' which given master fields 
for the functional integral with action 5'i + 5*2 + . . . allows us to construct master fields 
for the combined functional integral. The multi-matrix (|2.3| ) is simply n copies of a single 
Gaussian master field acting on separate parts of the free Fock space. To generalize this 
we need more general one-matrix master fields. 

An explicit master field for a general one-matrix integral 

^^g-w(M) (2.7) 




* 



with a similar form is given in 



M 



«+E 



i2.i 



Since the coefficient Cm only affects the expectations Tr M"^ for n > m, clearly we can 
find Crn's which reproduce any W. The diagrammatic argument works the same way: we 
associate a diagram to each term in the expansion of (|2.2|) , with each power of a* creating 
an additional internal leg. The free nature of the operators prevents lines of this diagram 
from crossing and thus the diagrams are planar. 

The coefficients Cn are associated with connected diagrams with n + 1 legs in this 
argument and thus we can write 

M = a + f{a*) (2.9) 



where / is the generating functional of connected planar diagrams 



/H = E^Si 



w 



(2.10) 



n>l 



as introduced in Brezin et. al. Although many generalizations of ( p.8|) to more general 
matrix models can be made, this observation (due to R. Gopakumar and D. Gross |]1( 



will motivate the particular one we introduce below. Here, it means that f{w) is given by 
the results of [Q : they define the generating functions of disconnected diagrams 

(PU) = V w"M^(m") = W( — i — ) (2.11) 

■^^ 1 — um 

n>0 

and connected diagrams 

v^(j) = i + $^Gi^)r 

n>l (2.12) 

<Piu)=ijU) (2.13) 

j = ui^ij). (2.14) 



The relation is then 



with j defined implicitly by 



In words, every time we attach to an external leg with u, we allow for the possibility of an 
arbitrary 'rainbow diagram' to its left. The sum of these is generated by expanding the 
implicit relation. 

For many purposes it is more convenient to work with the resolvent 

R(z) = W(^^) = V Vr(m")2-i-" = z-^iz-^) (2.15) 

z — m ^-^ 

n>0 

and the symbol of the operator M 

Kij)=r' + fij)=r'i^{j) (2.16) 

The relations ( |2.11| )-( p.l4| ) are then equivalent to |]I|] 

K{R{z)) = z. (2.17) 

The free product of several such operators, giving the master field for a decoupled set 
of matrix integrals, can now be constructed simply by taking [^ 

M^=a^ + Y,icc.)n{a*J'' (2.18) 

n 

with the operators as above. 

A virtue of this representation for present purposes is that the operation '+' between 
two free master fields (the 'additive free convolution' of Voiculescu) is very simple. Define 
a master field by 

{M^+2) = (Tr (Ml + M2)"). (2.19) 
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If Ml and M2 are relatively free - i.e. they act on independent factors in a free product 
of Hilbert spaces - this is completely determined by Mi and M2 • We can replace 

Ml +M2 = ai + a2 + /i(a*,a*) + /2(a*,6*) (2.20) 

by 

Mi+2 = ai + /i(at,a*) + /2(at,6*) (2.21) 



because in each term in the expansion of ( p.2|) , only a single one of ai or 02 could contribute. 
In diagrammatic language, an appearance of Mi + M2 produces a sum of terms where we 
attach a diagram from Mi or M2, and we do not have to keep track of which of the two 
the diagram came from. The simplest statement of this for a single matrix problem is 
in terms of the function f{j) of ( p.9|) : 



fi+2iJ) = fiU) + f2U). (2.22) 

In stochastic quantization, the noise t^ir) at each time r is free with respect to every- 
thing at previous times, and thus this result will apply to the Langevin equation for the 
master fields. 

A natural generalization to multi-matrix models is to consider representations of the 
form 

Mc.=a^ + Ual,a*2,...,a*J (2.23) 

which we will call a master field in the 'C-representation.' The functions /q, are general 
functions depending only on the creation operators. Since these operators do not commute, 
these are functions on a free algebra and a basis for these is rather large: for example, 

/(a*, b*) = Ma*) + /g(6*) + fabianf^iibn + /fe„(6*)/,a(a*) + . . . . (2.24) 

As for the one- matrix c^'s, there are enough coefficients in a series expansion of these 
functions to reproduce functionals Ws for very general actions. (The question of which 
Ws can be reproduced is non-perturbative.) 

From the matrix integral, the result ( |2T^ ) would be computed as a sum of disconnected 
diagrams and thus, following the observation of Gopakumar and Gross, we reproduce it 
if fa is the generating function for connected diagrams with a marked external Mq, leg, 
with a term a* ... a* creating internal legs M^^ . . . M^^ . In writing the expression ( |2.2[ ) we 
break the cyclic symmetry of the trace in the original Green's function. This allows us to 
assign each connected component of a disconnected diagram a unique 'first leg,' the first 



to appear in reading (|2.2| ) from right to left, and thus there are no additional symmetry 
factors: each disconnected diagram is generated once in the expansion. 

The result of Brezin et. al. generalizes to any number of matrices. We again introduce 
generating functions ^(wa) and i/j{ja) of disconnected and connected diagrams respectively. 
The variables u^ and j^ are free (non-commuting) and we have expansions like 

<P(u) = 1 + ^ <M^(m^) + Yl «2 W^(^2 ) + --- + Y1 <u^W{m'['m^) + ... 

1 (2.25) 

= 'W( ^^)". 

1 — um 



The analog of (gT3| ) - ([2Tl is 



(2.26) 



These sum all rainbow diagrams. 

The Gaussian two-matrix case is a good illustration: 'il}{j,k) = 1 + p + k^ so j = 
y + y{f + k^), k = z + z{f + P) and ^(y, z) = I + y'^ + z"^ + y{f + P)y + y'^{p + /c^) + 
y(i^ + k'^)yij'^ + ^^) + ■ ■ •■ The grading of the terms by word length makes this expansion 
well defined. 

Can this relation be inverted to compute i/^ij) from (p{u)7 Although there is a theory 
of functions of free variables (e.g. W^), we know of no general treatment of such problems. 



We want to interpret (p.26| ) the other way around, solving j^ = Uct4>{u) for w as a function 
of J. Writing z^ = w^^ and Za = 4>{^~^)Ja^ makes the inverse relation very similar to the 
original, as with (|2.17| ). 



Whether we can solve for ij: given cf) in practice, whether there is any analog of the 
analyticity properties of R{z) and K{z) in the complex plane, and so forth, is unknown. 
Making less formal definitions of these generating functions would also be valuable to 
justify the claim (which we believe - of course this is well understood in the one matrix 
case) that the objects under discussion are well-defined beyond perturbation theory and 
that all references to diagrams in the preceding are pictorial only. 

The work of [jl| does not assume that the variables are hermitian matrices and one 
can treat gauge theories in this formalism either directly in terms of the A^ as above, or 
by taking holonomy variables U = expiA as free variables with spectral density on the 
unit circle in the complex plane. The latter would be appropriate for lattice gauge theory 
or the principal chiral model. 

It is worth pointing out a simple application of freeness (in Voiculescu's sense) in 
this context. In the extreme strong limit of lattice gauge theory, the link integrals are 
completely independent, and in the large A^ limit they go over to free random variables. 



The master field is just Ui -^ Ui with each Ui a generator of a free group and tr O = 5{0,1). 
One expects that this is the way confinement will be realized in any large A^ gauge theory:^ 
gauge-invariant correlators of non-trivial operators such as tr F{x)UxyF{ii)Uyx will fall off 
exponentially with distance and F[x) and F{y) will become relatively free. By working in 
axial gauge {x—y)-A = and integrating, it is fairly easy to argue that {A{x)A{y)) ~ \x—y\ 



and that this will give confinement. \lt 



3. Stochastic Quantization 



This is well treated in the literature []T2| and we just mention a few of the main points. 
The Langevin equation is 

--(/)a(r,x) = -— — -S[(i)\+r]^{T,x). (3.1) 

OT d(pa{T,X) 

A quantum expectation is reproduced in terms of the large r limit of a solution 0a (t) — * (pa 

as 

(<?^ai<^a2 ••■) = (<^ai<^a2 •■•)r; (3.2) 

where the new average is over Gaussian noise with 

{Va{r)vp{r')) = 26apS{r-r'). (3.3) 

The probability distribution P[4>, r] satisfies the Fokker-Planck equation 

_^P.5^ * («p)+ilp (3.4) 

dr ^ S(/}a \S(t>a J m 

with large r limit -P[0] = exp — S'[0]. 

Perturbatively one can describe the result as a sum of 'stochastic diagrams.' The 
stochastic diagrams for (0i . . . (pn) are formed by taking all tree diagrams contributing to 
each cj)ct (one external leg) with sources rj (zero or more external legs), and tying these 
together in all possible ways with a quadratic rf' vertex. The Gaussian case illustrates the 
role of the 2 in ( |3.3| ): the single diagram has two propagators e~'^^^ "'""^ -* and a vertex. We 
integrate the result j dr, and the vertex factor 2 compensates the extra 1/2 produced by 
the integral. 

Recently a nice interpretion of the proper time in these diagrams has been made by 



Jevicki and Rodrigues and by Ishibashi et. al. [|0|- it is the analog of the time coordinate 



This argument was developed in a conversation with D. Minic. 



r on a string world-sheet in 'temporal gauge.' In this gauge, r is defined as the minimum 
geodesic distance to a boundary of the diagram. 

All loop diagrams come from sewing tree diagrams with rf vertices and thus we can 
regulate the theory by regulating the {rjrf) correlator. This is 'stochastic regularization' and 
is particularly attractive for gauge theories. [1l4| , p!5[| The present discussion will however 



assume that the theory of interest is already regulated, for example by the lattice. 

4. Stochastic Master Fields 

Greensite and Halpern pointed out that when the variables ?7q-(t) are large A^ matrices, 
a 'master noise' can be constructed to reproduce (|3.3|) . In the present language this is just 



I) . This leaves us with the problem of defining the rest of (|3.1| ) . 
We can define the first order proper time derivative in ( p.l|) as the e — > limit of the 
difference equation 

M^{t + e) = M^{t) + eM^{T). (4.1) 

We will assemble MQ,(r) as the sum of two terms. One term is ?7(t), which is always free 
with respect to the rest of the equation. We can thus rewrite the equation 

M = riir) (4.2) 

using {^^) as 

Ua*) = 2al. (4.3) 

In the one-matrix case this is equivalent to 

kij) = 2j. (4.4) 

Using dr{K{R{z))) = 0, we can change variables in this to R{z), and find p,0 

R{z) + d,R{zf = 0, (4.5) 

the Hopf equation of matrix quantum mechanics. Its general solution is thus most simply 
written in terms of K{j). 

The forcing term in the equation is harder to write. Let us first do it for the one- 
matrix case. It is most simply written in terms of R{z) because we can derive it by varying 

M: 

Tr ^ — - = Tr ^-—rSM(M) 

z-M -5M{M) [z-Uy ^ ' 

= -a.[Tr^5M(.)]_. 
10 



where [-F]- is the projection on the part F{z) ~ 1/z at infinity, i.e. 

with the contour around the cut in F{z'). Thus we can rewrite 

M = -V'{M) (4.8) 

as 

R{z)=d,[V'iz)R{z)]_ (4.9) 



The imaginary part p{z) = {l/n)luiR{z — ie) is the spectral density of M and thus ( ^^ ) 
is essentially the analytic continuation of the first term in ( p.4|) . 

The full stochastic evolution is just the sum of the two terms. We can change variables 
to R{z) and find 

R{z) = d, [V'{z)R{z)]_ - d,R{zf. (4.10) 

As r ^ cxD, the solution approaches a stationary master field with R{z) = and (integrat- 
ing once) 

[V'{z)R{z)]_=R{zf + c. (4.11) 

A solution with R{z) ~ 1/z exists only if c = and this is the usual factorized Schwinger- 
Dyson equation for the one-matrix model. 

The new possibility is to change variables in (|4.9| ) to K{j) to find 



K{j) = -dj \j V'{K{j))]_ + 2j (4.12) 

and at the stationary point, 

dj [j V'{K{j))]_ = 2j (4.13) 

which we can integrate to get 

j-^[jV'{K{j))]_=j. (4.14) 

We can regard it as a precise version of 

V'{M) = ry, (4.15) 

a classical equation of motion for the master field. 

Care must be taken in interpreting (|4.15| ) as a naive treatment of it, for example using 
it to conclude that we can derive the spectral density p(A) at the saddle point by the 
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simple change of variables V'{X) = x from a semicircular p(x), is incorrect. We see that 
the calculus of [|T| provides a correct interpretation. 

It is interesting to ask whether we can interpret the equation ([4.15| ) more directly 



as an equation for the operator M. Strictly speaking, at the fixed point r -^ oo, we 
have not M = but rather that M is an infinitesimal similarity transformation ('gauge 
transformation') of M. We should write 

= M = -V'{M)+r]+[A,M] (4.16) 

for an operator A to be determined. 

If we compute V'{M) = ^ kVkM^~^ in the C-representation, by defining M'^~^ as a 
power of ( |2.8D , the result will not be of the form a + /(a*). We need to choose A to fix 
this. 

For V'{M) = m?M, the time evolution M = -V'{M) is 

M + eM = {1- em^)a + (1 - em^) ^ Ck{a*)''. (4.17) 

k 

This is not in the C-representation but we can do a simple 'canonical transformation' 

a ^ (1 + em'^)a and a* — *> (1 — em^) a to find 

M-[A,M] = -m2^(A; + l) Cfc (a*)^ (4.18) 

k 

This is precisely the time evolution of ( ^.12| ). We see that both —V'{M) and 77 must be 
regarded as infinitesimal variations, and in this sense ( [4.16| ) is valid. 

More generally, we need a similarity transformation which eliminates higher powers of 
a. This is a bit more complicated, as is suggested by the operation [j V'{K{j))]- in ( ^4.14| ) 
which is defined to eliminate positive powers of u, not j. Let A r^ B indicate that two 
operators differ only in matrix elements Aij — Bij with i and j less than some finite integer. 
Then, [M,a] ~ [M,a*] ~ 0. We need [[A,M],a] ~ and this implies [A, a] ~ p{M) for 
some function p. This will produce [A, M] ~ p'{M) and the result (which can be verified 
by comparison with ( [1.11|) ) is that ( [4.16| ) becomes 

V'iK)--p'{K)=j (4.19) 

J 

where p'{K) is the unique polynomial in K which cancels negative powers of j from V'{K). 
For a multi-matrix model, we have 



m = 2j2jl (4.20) 
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The forcing term is again simple in terms of disconnected diagrams: 

a,(Tr M^^M^^ ...) = - 5^(Tr M„, . . . d^^S[M] ...). (4.21) 

n 

The derivative 9qMq,^ . . . here is a cychc sum of terms 5a,aiMonj^^ . . . Mq,^_^. The diagram- 
matic interpretation is a sum of terms where we attach k — 1 legs of S[M] to an existing 
diagram, leaving one leg. 

If we change variables in (|4.20|) and go to the fixed point, this will produce the standard 



'integrated' Schwinger-Dyson equation: 



"y/J^^^T^ M^, ...M„ ) =y"( ,-%TrM«, ...M« .) (4.22) 



In a sense this is the multi-matrix generalization of the change of variables from ( |4.4| ) to 
( ^75| ) which solves the Hopf equation. In computing this, the relation j^ = Uai^ gives us 
if) = 2 ^(mq,(^)^ which is the piece of the right hand side of ( |4.22| ) in which a derivative acts 
on the first matrix; the rest comes from writing (/>(w) = V'(i) but taking into account the 
variation of i/j in the relation j^ = Uaifj- This must produce the other insertions of u^ ■ ■ ■ Ua 
necessary to produce a cyclically symmetric answer, and we verified this diagrammatically. 



5. Comparison with Previous Work 

After the completion of the original version, we discovered the almost forgotten papers 



m and 18 



In a beautiful and prescient work of Haan [|T8[ , the first precise master fields for higher 
dimensional large N matrix models were constructed. An explicit construction for the D- 
dimensional Gaussian field theory essentially identical to ( |2.3| ) and (|2.4|) was given. It 



was shown for the field theory with Tr M"* interaction (by an argument which generalizes 
trivially to any action) that the master field satisfies the equation 

.^^+2n.|) = (5.1) 

dMi 

with 

[U,,Mj] = -id,j |)(|. (5.2) 

(This algebra arose independently in @,|[). The equation was studied explicitly in D = 
and for Gaussian field theory. (See also fl^]-) It was argued that the usual adaptations 
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of the GNS construction for quantum field theory demonstrate that such a master field 
exists. 

In [jl^, Cvitanovic, Lauwers and Scharbach gave the relation ( |2.26| ) between the gen- 



erating functional of planar disconnected and connected diagrams Z'[j] and W[J] (treated 
as functions of free variables), as well as the relation with the generator of IPI diagrams. 
A number of useful results were derived: one was the following representation of the dis- 
connected planar Green's functions (their equation (4.5)): 



fdW[J] d 



^----H'w-wyi^^ 



dJ, 



1, (5.3) 



where the J^'s are free variables and 

_d_ 

If we make the identification 



Jjf{J)=d,jf{J). (5.4) 



a* = Ji] ai = — (5.5) 

we see that this is identical to ( p.23|) . They also derived the equation of motion 

1 = Ji- (5.6) 



' S 



dM, 



dW[J] d_ 



dJ dJ 



However, they did not realize that these operators could be interpreted as a master field. 



6. Conclusions 

The theory of operator algebras allows one to construct master fields for general large 
N field theories. A particularly useful version of this for special cases can be taken from the 
work of Voiculescu. Gopakumar and Gross have recently explained its relation to conven- 
tional diagrammatic methods: this master field encodes the generating functions for con- 
nected planar diagrams. Using this observation it is straightforward to adapt Voiculescu's 
construction to general large N field theories. We studied the stochastic equation of Green- 
site and Halpern for the master field, and explained the relation of the equation in this 
representation to the standard factorized Schwinger-Dyson equations. 

Perhaps the most important outcome is that a number of seemingly disparate ap- 
proaches to large A^ have now been connected in a precise way, and the subject put on 
a firmer footing. Along the way we rediscovered the neglected but we believe important 



work of [0 and [|T^ . In particular, (|5.6|) appears to be the simplest form for the equation 



of motion determining the master field. 
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The connection with the work of Voiculescu is recent and many of his results and 
derivations are new for the large A^ application (such as the derivation of matrix quantum 



mechanics of H; see [O] for more examples.) His concept of freeness is fundamental and 



can be used in a number of physics applications - an example is the Neu-Speicher model of 



localization. |21|] In section 2 we pointed out its relevance for confinement. More generally, 
the relevance of the theory of operator algebras to large N is now clear. 

The concept of master field is rather protean and we believe it will be useful to consider 
other representations (equivalent up to similarity transformations) as well, and be able to 
convert between representations. We might draw an analogy with the use of different 
gauges in analyzing gauge theory to exhibit different properties of the theory. The GNS 
construction directly relates the original Green's functions (|2.1D to another master field; 



thus a simple example of a problem in this direction would be to show that the relation 
( p.26|) is invertible in practice, determining the master field in the C-representation given 
the original Green's functions. Probably, other valuable representations for multi-matrix 
and higher dimensional models remain to be discovered. 

To test these ideas we need to work with higher dimensional theories. Even matrix 
quantum mechanics (with a general potential) qualifies, since an explicit master field de- 
scribes non-equal-time correlators, and thus contains more information than the existing 
solution. 

At present, we would guess that the question of whether or not a particular higher 
dimensional theory is integrable or not has the same answer at finite A^ and large A^. If 
so, applications to non-integrable theories such as higher dimensional gauge theory will 
require developing approximate methods. 

Developing analytic approximations for loop functionals and master fields is an im- 
portant open problem, and it is here that we believe the theory of operator algebras will 
show its deeper relevance for these problems. Paradoxically, we did not really use it in 



the present work: the algebra defined by (|2]J) was treated in a purely formal and com- 



binatoric way. In reality the theory of operator algebras is a branch of analysis, and our 
hope is that out of it and the physical applications will eventually grow practical tools 
which allow representing realistic loop functionals and master fields with a finite amount 
of information. 

We thank R. Gopakumar, D. Gross, M. Halpern, M. Li, D. Minic, G. Moore, S. 
Shenker and I. Singer for useful discussions. We especially thank M. Halpern for pointing 
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After the completion of the original version we received the work of R. Gopakumar 
and D. Gross [0, where the representation (|2.23| ) is also proposed. 

This research was supported in part by DOE grant DE-FG05-90ER40559, NSF PHY- 
9157016 and the A. J. Sloan Foundation. 
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